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PERFORMANCE AND STABILITY ANALYSIS OF 


A PHOTOVOLTAIC POWER SYSTEM 
by Walter C. Merrill, RonaIdJ. Blaha, and Roy L. Pickrell 
Lewis Research Center 

SUMMARY 

Recently much emphasis has been placed on the conversion of solar energy to more 
useable energy forms. One method is the direct conversion of solar insolation to a dc 
voltage and current by the use of solar cells. This dc power can then be "inverted" to 
ac power compatible with commercial utility lines to drive various ac loads. This 
report studies the performance and stability characteristics of one such system, a 
10 kVA photovoltaic power system. The system is studied using linear Bode analysis 
and a nonlinear analog simulation. Power conversion efficiencies, system stability, 
and system transient performance results are given for system operation at various 
levels of solar insolation. Additionally, system operation and the modeling of system 
components for the purpose of computer simulation are described in this report. 


INTRODUCTION 

Recently much emphasis has been placed on the conversion of solar energy to more 
useable energy forms. One method is the direct conversion of solar insolation to a dc 
voltage and current by the use of solar cells. This dc power can then be "inverted" to 
ac power compatible with commercial utility lines to drive various ac loads. This 
report studies the performance and stability characteristics of one such system, a 
10 kVA photovoltaic power system. 

The Terrestial Photovoltaic Project Office, as part of an DOE effort, is developing 
a 10 kVA Inverter/Controller subsystem for testing in a Photovoltaic Power System 
(PPS), at Lewis Research Center (LeRC). The purpose of the combined system (PPS) 
is to investigate the application potential of photovoltaic power sources to supplement 
and/or replace conventional energy supplies. The PPS is designed for operation in 
either of two modes, the stand alone mode (mode I) and the utility mode (mode II). 

Mode I requires a system load independent of other power generation sources. Power 


injection in this mode is a function of load only. In mode II power from the PPS is in- 
jected into a commercial utility distribution system. The amount of power transferred 
in this mode will depend upon the solar array and its environmental conditions. Mode I 
operation is not considered in this report. 

The PPS consists of five major blocks: (1) an array of solar cells, (2) the 10 kVA 
dc to ac power inverter, (3) the power system controller (PSC), (4) the pilot sensor 
array, and (5) the utility interface . The solar array directly converts solar insolation 
into dc power. The inverter transforms this dc power to 60 Hz ac power compatible 
with most electrically driven machines. The inverter design is an adaption, under 
LeRC contract, of an existing product line. Included in the inverter design is the in- 
verter phase control, which adjusts the inverter output voltage phase to assure that the 
power injection into the utility line is at unity power factor (i.e. , the output current is 
in phase with the utility line voltage). The PSC monitors system operation and adjusts 
array operating voltage (by setting the power level of the inverter) to achieve maximum 
power operation of the array. The power system controller (PSC) conceptual design 
was accomplished at LeRC. The detailed design, fabrication, and acceptance testing 
was accomplished under LeRC contract. The pilot sensor array provides the reference 
voltage from which the target operating voltage for the solar array is derived. The 
utility interface supplies a power sink, for absorbing the PPS output power without dis- 
turbing the utility voltage magnitude or phase (in operations and analysis the utility 
voltage is used as the reference phasor) . 

This report begins with a description of the major elements of the PPS. The re- 
port then describes an analysis of the stability and overall performance of the PPS. In 
particular, a linear stability analysis of the system at a high and low power operating 
condition is discussed. Additionally, a nonlinear analog simulation to more completely 
model component performance and system stability is presented. Finally, a discussion 
of the analysis results concludes this report. 

PPS DESCRIPTION 

For the subject simulation and analysis, the PPS is defined in block diagram form 
in figure 1 . Dynamically, the important elements are (1) the solar array, (2) the power 
system controller, (3) the power inverter with its input and output filters (a 0 . 021 
Farad capacitor and a 0.005 Henry choke, respectively), (4) the inverter phase con- 
troller, and (5) the reference utility line voltage. In this equivalent circuit, the pilot 
sensor array function is obtained from the solar array open circuit voltage data model. 
Figure 1 also indicates important system variables and the information/command flow 
in the system. A list of symbols is given in ajpendix A. 

Energy conversion is accomplished by the solar array. Impinging direct solar 
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insolation is changed to a dc voltage and current. The power inverter converts dc 
power to ac power. The input filter capacitor isolates the array from the pulsating dc 
input current requirements of the inverter. The inverter output filter inductor sup- 
presses the high frequency harmonics present in the synthesized 60 Hz output voltage. 
The PSC (in mode II) sets the operating voltage of the solar array by controlling the 
magnitude of the inverter output (generated) voltage. The inverter phase control main- 
tains an inphase relationship between output line current and utility line voltage by ad- 
justing the phase of inverter generated voltage. Each element is now discussed in 
more detail. 


Solar Array 

The solar array consists of several himdred interconnected solar cells. The 
current-voltage characteristics for each cell for different solar insolation levels and 
ambient temperatures are determined experimentally and are presented in figure 2. 

The entire array consists of 64 parallel connected strings of solar cells with 600 series 
connected cells per string. Therefore a total solar array characteristic could be ob- 
tained by a simple scaling of the individual cell characteristic as shown in figure 2 . 

The ordinate intercept value of each curve is called the short circuit current and the 
abscissa intercept value the open* circuit voltage. A normalization of each curve by its 
open circuit voltage and short circuit current values yields a single curve representa- 
tive of all temperature and insolation conditions (fig. 3). This representation is con- 
venient for usage in a system simulation since only one nonlinear curve need be pro- 
grammed. 

The solar array characteristic of figures 2 and 3 can be approximated closely by 
a diode equation (ref. 1). 



The constants were graphically determined for the data presented as A = 10 and 
B = 1. Solar array power is then 


^SA “ ^D%A 


(2) 


The maximum power point is found from 


3 


( 3 ) 


3P, 


SA 


= 0 


3E 


D 


From (1), (2), and (3) the necessary optimality condition is 

Ep _ exp 10(1 - Ep/Ep^ ~ ^ 
^D, N 

This equation is satisfied when 


( 4 ) 


Ej^« 0.7823 (5) 

Thus, when operating in the presence of direct insolation, the array generates a dc 
power in accordance with the voltage-current characteristics of figure 2. Maximum 
power operation at all conditions can be approximated if the array output voltage is 
maintained at 0.80 times the open circuit voltage. In mode II operation this is one 
goal of the PSC . 


Power System Control 

The mode II configuration of the PSC is given in block diagram form in figure 4 . 
The feedback voltage, Ej^, is filtered by a three pole, 40 Hz, Butterworth filter to 
remove 120 Hz noise and compared to a reference voltage, Ej^p. The voltage, 
^REF’ is selected proportional to Eq^^ such that 

Ep ^ 0.8 Eq^ 


when 


^D,F ■ ^REF " ® 

The error signal is filtered by a 500 Hz, single pole, low pass filter to remove high 
frequency noise. The filtered error signal then drives a proportional-integral control 
to generate the inverter control voltage . 

Power Inverter 

The power inverter transforms dc input power generated by the solar array into 
single phase ac power compatible with most electrical machines. Conversion is 
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accomplished by approximating the sinusoidal inverter output voltage by a prepro- 
grammed sequence of dc pulses. The resultant sequence is filtered to remove high 
order harmonics. Conceptually, the inverter can be considered an ideal power con- 
version device, that is, dc power in equals ac power out at every instant of time. Ac- 
tual inverter losses can then be approximated by two resistances on the output side of 
the inverter. Figure 5 is a schematic of the filter capacitor, power inverter, and in- 
verter output impedance. The resistance loss of the filter choke is negligible- Fig- 
ure 6 is a phasor voltage diagram that approximates the dynamic magnitude and phase 
relationships of variables given in figure 5. Note that a steady state condition would be 
represented by a right triangle diagram where ^ = 0. Also, phase is defined with re- 
spect to the utility line reference voltage. 

Inverter control is accomplished by two means. First, adjusts the mag- 

nitudes of the generated voltage and the line current. These magnitudes also are re- 
lated by the inverter output impedance. The control variable, 6, sets the phase angle 
of the generated voltage to control the phase difference between line voltage and cur- 
rent. This phase angle is adjusted by the inverter phase control (IPC) circuit. 

Inverter Phase Control 

The inverter phase control, IPC, determines the phase difference between line 
voltage and line current by comparing square waves derived from these two signals. 
This difference serves as an error signal to a proportional-integral control. The con- 
trol output drives a voltage controlled oscillator which can affect minor adjustments in 
the frequency of the generated voltage (ref. 2). These adjustments effectively set the 
generated phase, 6 . The IPC, therefore, acts as a phase locked loop that sets invert- 
er generated phase such that (p = 0 in steady state. This inphase line voltage-line 
current condition insures maximum power injection by the inverter into the utility 
power line. 


LINEAR ANALYSIS 

In the initial phase of the investigations, a linear analysis was undertaken to in- 
vestigate the stability of the PPS. This analysis was required to verify basic stability 
of the total system since individual components were being manufactured at different 
locations. Several assumptions were made to simplify this "first look" analysis. 
First, the system was linearized at both a maximum and minimum power operating 
condition to obtain Laplace transfer functions for each of the blocks in figure 7 . Sec- 
ond, the inverter was assumed to be in a phase-locked mode at all times so that only 
the power control loop was active. Standard Bode analysis (ref. 3) was performed 
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using the transfer functions of appendix B. The dynamic resistance of the solar array, 
Rg^, is defined as the negative reciprocal of the slope of the solar array characteris- 
tic at the point of maximum power. For a high and low level of insolation two values 
of Rg^ were used, 

RgA(high) = 4 
RgA(low) = 200 

The results of the Bode analysis, given in figures 8 and 9 demonstrate that the 
nominal or high insolation condition is only marginal stable . A stability condition 
exists when at 180° phase shift the open loop gain is less than 1. At the low insolation 
condition the system is unstable by a gain margin of 8 . 5 dB (gain factor of 2 . 7 ) . A re- 
duction in closed loop static gain of 12 dB will stabilize the system at both operating 
conditions leaving adequate stability margin at the low insolation condition. Thus, the 
PPS configuration with the reduction in loop gain can be expected to operate in a stable 
manner over a wide range of insolation levels. Further analysis of the impact of sys- 
tem nonlinearities and the IPC on system performance and stability was desired. To 
complete this analysis a nonlinear, analog simulation of the PPS was devised. 


NONLINEAR ANALYSIS 

The nonlinear analysis of the PPS further evaluated the stability of the system to 
varying solar insolation conditions. Also, transient system performance and tolerance 
to line voltage disturbances were studied. This section describes the nonlinear simu- 
lation and then discusses simulation test results. 


Nonlinear Simulation 

The nonlinear simulation was devised by combining component models to simulate 
the information flow of figure 1. Equations for each component are given in appen- 
dix C . The normalized solar array characteristic was approximated by 15 straight 
line segments to simulate solar array operation. Four characteristics generated by 
the analog simulation from the normalized curve by impropriate scaling are presented 
in figure 10 for comparison with the given characteristics of figure 2 . Modeling equa- 
tions for the PSC and the filter capacitor were determined from their circuit diagrams 
in a straightforward manner. Conceptually, the inverter was modeled as in refer- 
ence 4. Equations and parameter values for the inverter and inverter phase controller 
components were obtained from the inverter manufacturer. Equations that model or 
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approximate the dynamic operation of the inverter output filter impedance were deter- 
mined from trigonometric relationships associated with the phasor diagram of figure 6 . 

Equations (C3) to (C5) are a macroscopic or input-output representation of invert- 
er operation. Considerable detail about internal inverter operation, such as the actual 
dc to ac voltage inversion and control processes was eliminated to simplify the simula- 
tion. Also, some approximations were incorporated such as considering the phase of 
Eji^ and E^ to be identical. Digital computer simulation of the exact trigonometric 
equations at several operating points showed the approximation error to be negligible. 
In equation (C14) the approximations 


<p = sin cp 


9 = sin 6 


are used. Equation (C14) should be written with and 9 as dependent and indepen- 
dent variables, respectively. However, since cp and 9 are generally less in magni- 
tude than 0.7854 radians (45°), and since sin 9 and not 9 is available from the simu- 
lation, the approximation of (C14) was used. 

Another simplifying approximation is the phasor triangle representation of figure 6 
for inverter output impedance. This is a more important limitation since it limits the 
maximum useful simulation bandwidth. The voltages and currents in the impedance 
network are usually sinusoidal in nature. In this case the phasor representation of fig- 
ure 6 is well defined and accurately models the process. However, transient results 
for time intervals less than one period (0.0159 sec) of the phasor frequency (60 Hz) and 
disturbances that generate a nonsinusoidal inverter output voltage are not accurately 
represented. These modes of operation can occur for brief intervals during line tran- 
sients. Because of their brief duration and because the simulation accurately models 
the trends for these brief intervals, this limitation does not seriously impact the use- 
fulness of the simulation. 


Experimental Results 

Several experiments completed with the nonlinear simulation were used to evaluate 
system (1) stability, (2) transient performance, and (3) steady-state efficiencies. 

Stability 

The linear analysis indicated that instability, if it were to occur, would be most 
likely at low insolation levels. Thus, to determine the critical insolation level at 
which instability occurs, as well as simulating system operation during a worst case 
"cloud cover" transient, the following experiment was devised. Insolation level was 


7 



o o 

decreased at constant rate from 100 to 6.5 mW/cm in 5 seconds, held at 6.5 mW/cm 
for approximately 6 seconds, then increased at constant rate to 100 mW/cm in 5 sec- 
onds. The experiment was conducted using the nonlinear simulation with a constant 
line voltage magnitude, = 340 volts, and a constant reference point setting of 

Ejj = 0.8 Eqq- The PSC forward gain was nominally selected as K = 5. This settii^ 
corresponds to the reduced stabilizing gain found in the linear analysis. Insolation 
level changes were simulated by varying the nonnalizing values for the solar array 
characteristic according to (C19). 

Time records of the cloud cover transient are given in figure 11 for several sys- 
tem variables. For the cloud cover transient ambient solar cell temperature was as- 
sumed to be 55° C. Overall, system operation is stable for almost the entire tran- 
sient, that is, a nonoscillatory, smooth trajectory for the solar array operating point 

is achieved. However, a limit cycle is observed in the inverter phase control portion 

2 

of the system for insolation levels less than 20 mW/cm . Additionally, the solar 
array voltage, Ej^, oscillates about its commanded value. This loss of control and 
subsequent inverter phase control oscillation is a result of the small difference in 
magnitude between generated and line voltages. From figure 11 the onset of the limit 
cycle occurs when E^ jyj (or j^) appears equal to E j^. From figures 5 and 6 
Ef M always must be greater than Ej^ in steady state. From table I at 
20 mW/cm^, Ej^ jyj - Ej^ = 0.72 volts. This is a difference of only 0.18 percent. 
Practically, the analog computer used to implement the nonlinear simulation cannot 
resolve this difference . Conceptually, this can be modeled as a deadzone nonlinearity 
in the system. Therefore, ability to set 6 and injected power is lost, a power im- 
balance is incurred, and localized limit cycling results. Similar results may be anti- 
cipated from the actual system since it will be equally difficult for any real hardware 
to maintain accurate power level control Ity resolving such small voltage differences. 

Transient Performance 

The transient response of the PPS to changes in the commanded reference point 
and to changes in the utility line voltage were also simulated. Responses to reference 
point step changes define dynamic response times of the system. System response to 
utility line voltage fluctuations is particularly important since this is typically the most 
severe operating constraint on utility mode operation. 

Figures 12 to 15 show system response to ±10 percent changes in Ej^p for four 
insolation levels, 100, 33, 14, and 6.5 mW/cm^, respectively. These experiments 
were conducted with Ej^ = 340 volts and K = 5. As Ej^ chaises the PSC and the 
inverter must adjust the generated voltage, line current, etc. to satisfy the change in 


8 


available solar power. Since Ej^ is constant, Ij^ must change in direct propor- 
tion to available power. From (C4) and are directly related in steady 

state. In each case the overall dynamic response, reflected chiefly in is 

the same. At the two low insolation level conditions the IPC cycling is evident. Al- 
though precise voltage control and power injection is not being accomplished in these 
two cases, global system stability is maintained in the presence of the operating point 
disturbances. In the 100 and 33 mW/cm cases adequate set point control is accom- 
plished. System response is basically characterized as imderdamped and second 
order. More damping, if desired, could be obtained by reducing PSC forward gain. 

Figures 16 to 19 show system response to ±10 percent changes in the magnitude of 

line voltage. These changes represent worst case line voltage transients that the PPS 

must tolerate. These experiments were also conducted at 100, 33, 14, and 6.5 mW/ 

2 

cm . Nominally, Ej^ = 340 volts and K = 5. As Ej^ step changes are encoun- 
tered, the IPC adjusts 6 and the PSC adjusts Eq^j^^ to match the new level of line 
current. For the 100 and 33 mW/cm^ experiments line transient tolerance is good. 
Notice that as insolation level decreases the line transients have a more pronounced 
effect on system operation . This is true since the 10 percent line voltage transient 
causes a higher percentage power imbalance at lower insolation levels. Again at the 
14 and 6.5 mW/cm levels, global system stability is maintained for these line tran- 
sients although accurate voltage and power control is lost. 


Efficiency 

This performance test was conducted to predict operating efficiencies and power 
levels at various solar insolation levels. Steady- state data were recorded and are 
listed in table I . System efficiency is calculated as 


1 ? = 


2EdIsA 



A comparison of simulation and design inverter efficiencies is given in figure 20. The 
values of the two power loss resistors in the inverter output impedance were selected 
to match the simulation to the design efficiencies at two power level points, no load and 
rated power. Resistor RLl represents resistive or load losses while RL2 represents 
no load power required for inverter operation. Figure 20 shows a good correspondence 
throughout the operating range. A comparison of the design and simulated power loss, 
as in figure 21, shows an essentially constant difference in power loss which the model 
does not account for. This loss discrepancy was not investigated in depth. However, 
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it was noted that in the inverter power switching transistors the diode junction losses 
are a linear function of load current rather than a squared function as indicated by the 
resistive element. Improved performance correlation could be obtained if necessary, 
by incorporating this source of power loss differently in the simulation. 


SUMMARY OF RESULTS 

A preliminary analysis of stability of the PPS demonstrated possible instability 
for the specified design value of PSC forward gain. From the linear analysis a gain 
reduction of 12 dB was found sufficient to stabilize the system at the high and low in- 
solation levels studied. Subsequent nonlinear simulation confirmed global system sta- 
bility over the entire range of insolation levels with the reduced PSC forward gain . 
However, a limit cycle was encountered in the IPC subsystem for insolation levels 
less than 20 mW/cm . This limit cycle phenomenon was found to be the result of a 
deadzone effect . The actual PPS will probably exhibit this limit cycle since the dead- 
zone is encountered for small differences in line and inverter output voltage which 
naturally result from low insolation levels. 

For insolation levels above 20 mW/cm the system exhibited (1) good operating 
point control of the desired solar array voltage and accurate control of utility line 
power injection for the cloud cover transient, (2) an underdamped, second order re- 
sponse to transients, and (3) acceptable rejection of ±10 percent line voltage 

step transients. For insolation levels below 20 mW/cm^ acceptable system stability 
and solar operating voltages were maintained for all the above experiments. However, 
oscillations in injected utility line px)wer and solar array voltage were encountered. 
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APPENDIX A 


SYMBOLS 

D( ) derivative operator, sec“^ 

E voltage, V 

exp Naperian exponential 

FI to F5 PPS transfer functions 

fl normalized array characteristic 

I current, A 

j 

K PSC forward loop gain 

L filter choke inductance, Henry 

R eiJective resistance, ohm 

RLl inverter loss resistance, ohm 

RL2 inverter loss resistance, ohm 

s Laplace variable, sec”^ 

7] system efficiency 

0 inverter output voltage phase angle, rad 

<p line current phase angle, rad 

Subscripts; 

A ideal inverter output (AC) 

CONT PSC 

D inverter input (DC) 

F filtered 

L line 

M magnitude 

N normalized 

OC open circuit 

REF solar array operating point 


SA solar array 

SC short circuit 

X filter choke reactive component 

1 actual inverter output (AC) 
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APPENDIX B 


F4(s) 


PPS TRANSFER FUNCTIONS 


Fl(s)- 


20 


0.000318 s + 1 


F2(s) = 


s + 1 


F3(s) = 


10 


0.159 s + 1 


0.0156 


(0.00794 s + 1) 


126/ \126 




F5(s) = 


^A 

0.021 Rg^s 


+ 1 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 
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APPENDIX C 


NONUNEAR SIMULATION EQUATIONS 


Solar array; 


Filter capacitor: 


^SA _ j 

^SA,N \ D,N, 




%A ~ 
0.021 


(Cl) 


(C2) 


Power inverter: 


^ _^A,M^A.M 


"D"D 


cos(0 - cp) 




^^^'CONT, F 
Inverter output impedance: 


j _ ^CONT ~ ^CONT.F 
0.159 


sin (p 


_ ^L.M ^X.M ~ ^1m 
M%,M 


(C3) 


(C4) 


(C5) 


(C6) 



(C7) 

^X,M 

(C8) 


cos 6 


= - si 


sin^e 


(C9) 
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cos (p - sgn(sin 6)\1 - sin <p 


^1,M “ ^A,M “ ^L,M 


(Cll) 


^A,M ^L,M 


(C12) 


cos(0 - (p) = cos 6 cos (p + sin 6 sin (p 


Inverter phase control: 


(D^ + 12.22 + 27.97 D)(sin 0) = -27.97(0.4368 D + l)(sin cp) 


Power system control: 


^REF ” 0-0125 


^D,N ~ ^OC 


ISA,N = 0.662 Ej^ j^- 164.05 


(CIS) 


(C14) 


(0.00794 D+1) (—\ +(^-^Vl (Ep» v-) = 0.0156 Et. (C15) 

L\ 126 / \ 126 / J ^ 

D jjo . 318(10 )D + ij “ K(D + 1)(E j, - E p) (C16) 


(C17) 


(CIS) 
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TABLE I. - SIMULATION PERFORMANCE DATA 


Variable 

2 

Insolation level, mW/cm 


65 

20 

38 

75 

100 

110 


Ambient temperature 




55 

55 

55 

55 

55 

0 

EOC* 

252 

266 

285 

292 

310 

375 

FsC’ ^ 

2.84 

12.0 

16.3 

32.0 

41.2 

46.4 

Ed. V 

202 

213 

228 

234 

248 

301 

ISA> ^ 

2.38 

10.4 

13.7 

26.8 

34.6 

38.9 

®L,M’ ^ 

340 

340 

340 

340 

340 

340 

^L.M- 

1.69 

11.7 

16.9 

34.5 

47.1 

63.8 

Powerg^, kW 

0.481 

2.22 

3.12 

6.27 

8.58 

11.71 

Power kW 

0.287 

1.99 

2.88 

5.87 

8.01 

10.85 

V 

0.596 

0.90 

0.921 

0.936 

0.933 

0.926 

e, deg 

0.536 

3.72 

5.35 

10.82 

14.6 

19.4 

®AM’ ^ 

340.96 

345 . 01 

347.50 

357 . 85 

366.95 

380.93 


340.01 

340.72 

341.49 

346.15 

351.35 

360.45 

^A,M’ ^ 

2.82 

12.87 

18.04 

35.60 

48.12 

64.63 

®CONT’ ^ 

0. 1120 

0.7791 

1.123 

2.289 

3.120 

4.216 


Power system control | 


3 pole 


\ 

Filter 

filter 


&PI 


Solar Isa Filter 
array ” cap 


Inverter phase control 


Figure 1. - Photovoltaic power system information flow diagram. 









































Figure 4. - Power system control block diagram. 
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e 

Figure 5. - Inverter schematic. 
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^t,M 

Figure 6. - Inverter output impedance voltage triangle. 



Figure 7. - Photovoltaic power system linearized block diagram. 
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Figure 12. - Concluded. 
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